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ABSTRACT 



We introduce a new quantity, the mass flux density of galaxies evolving 



from the blue sequence to the red sequence. We propose a simple technique 
for constraining this mass flux using the volume corrected number density in 



the extinction-corrected UV-optical color magnitude distribution, the stellar age 
indexes H5a and D n (4000) , and a simple prescription for spectral evolution 
using a quenched star formation history. We exploit the excellent separation of 
red and blue sequences in the NUV-r band Hess function. The final value we 
measure, p'r = 0.033 M Q yr -1 Mpc -3 , is strictly speaking an upper limit due to 
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the possible contributions of bursting, composite, and extincted galaxies. How- 
ever, it compares favorably with estimates of the average mass flux that we make 
based on the red lu minosity function evolution derived from the DEEPII and 
COMBO-17 surveys (bell et al.lbooibl : iFaber et alJliool ). p R = +0.034 M Q yr" 1 
Mpc~ 3 . We find that the blue sequence mass has remained roughly constant 
since z=l (p# ~ 0.01 M yr _1 Mpc -3 , but the average on-going star formation 
of psf — 0.037 M yr _1 Mpc -3 over < z < 1 is balanced by mass flux off the 
blue sequence. We explore the nature of the galaxies in the transition zone with 
particular attention to the frequency and impact of AGNs. The AGN fraction 
peaks in the transition zone. We find circumstantial, albeit weak evidence that 
the quench rates are higher in higher luminosity AGNs. 



Subject headings: galaxies: evolution — ultraviolet: galaxies 



Introduction 



There has been gro wing interest in the nature of the obser ved color bimodality in the 



distribution of gal axies (Balog h et a 



galaxy properties (IKauffmann et al. 



2004 : iBaldry et al.l 12004 ) , which is echoed in other 
20031 ). The color bimodality is revealed in a variety 



of color- magnitude plots, an d is particularly dramatic in the UV-optical color magnitude 
diagram (jWyder et al.l 120061 ) . The red and blue galaxy concentrations have been denoted 
the red sequence and the blue "cloud", although we elect to call both concentrations se- 
quences. Deep galaxy surveys are no w probing the evo lution of the r ed and blue seque nces. 
Recent work using the COMBO-17 (bell et al.ll2004bh and DEEP2 jFaber et alil2005h sur- 
veys provides evidence that the red sequence has grown in mass by a factor of three since 
z~l. It is natural to ask what processes have led to this growth, and in particular whether 
the red sequence has grown via gas rich mergers, gas-less (dry) mergers, or simple gas 
exhaustion, and whether AGN feedback has played a role in accelerating this evolution 
jSpringel. DiMatteo fc HernquislJliooi bi Matteo et al.ll2005h . 



We would like very much to identify galaxies which may be in the process of evolving 
from the blue to the red sequence. In this paper, we exploit the pronounced demarcation 
between blue and r ed sequences afforded by the Near UV-optical color-magnitude diagram 
Wyder et al.l (120061 ). We make the natural assumption that galaxies residing in the region 
between the blue and red sequences are making a transition from the former to the latter. 
We develop a metric (albeit model-dependent) for measuring the speed galaxies are evolving 
through this transition zone using spectroscopic indices. We use this to estimate the total 
mass flux of galaxies from the blue to the red sequence during the recent several gigayears. 
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Remarkably, the result we obtain is consistent with estimates based on the evolution of both 
the red and the blue sequences taken separately. 



We seek to determine the total rate that galaxies residing on the blue sequence evolve 
to the red sequence. Our approach is to make a series of simplifying assumptions, which we 
examine in section £j6j 

1. Galaxies evolving from the blue to red sequence pass relatively smoothly through inter- 
mediate colors. Galaxies do not "jump" , for example by hiding due to strong extinction 
only to be unvealed as a fully formed red sequence galaxy. In other words we assume 
that galaxies transition on timescales no shorter than the precision of the tracers we 
use below to determine those timescales. 

2. Galaxies can be described on- average in a given intermediate color-magnitude bin by 
a single parameterizable star formation history. In other words, we ignore composite 
stellar populations. 

3. We can ignore for the moment the detailed accounting of stellar mass growth due to 
merging. Rather, our goal is to determine the mass flux off the blue sequence, as a 
function of stellar mass, and then to total the mass flux leaving the blue sequence. 

4. The contribution of galaxies bluing (e.g., due to starbursts) from red to blue sequences 
through intermediate color bins can be ignored. 

With these assumptions, our methodology is quite simple. We define 0s(M r ) as the 
number flux per unit volume at a given absolute magnitude M r off the blue-sequence (unless 
explicitely stated the sign of all mass fluxes are defined to be positive). Given that we choose 
a narrow range of NUV-r color, there is a close correspondance between M r and stellar mass. 
Then the fundamental governing equation is 



where <f)(M r , NUV — r) is the volume-corrected NUV-r, M r color-magnitude function 
(Hess function), and r(M r , NUV — r; £) is the lifetime of the galaxy in this particular color- 
magnitude bin, for a star formation history £. 



2. Methodology 
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For most physically reasonable star formation histories, the lifetime in the bin is deter- 
mined by the rate of change of the color, such that: 

where y = NUV — r. dy/dt is positive due to assuption 4. The total mass flux is then 
h->R < Yl M ( M -) 4>(M r ,NUV-r) < ^[y;£] > (3) 

M r 

where M(M r ) is the stellar mass corresponding to M r . This equation is an equality only 
if all galaxies present in the intermediate color bin are moving from the blue to red sequence. 
Otherwise it represents an upper limit. 



We have determined <p(M r , NUV — r) in lWyder et al.l (120061 ). The key remaining issue 



is to determine dy/dt. Our approach will be to posit some simple quenched star formation 
histories and calculate stellar age indexes such as _D n (4000) for comparison with the observed 
galaxy indexes. These age measures will help constrain the family of possible quenching 
histories at the transitional colors. 



Data 



3.1. The Hess Function 



We generat e the y olume-corrected NUV-r vs. M r color-magnitude diagram as discussed 



m iWyder et al.l (120061 ) . We briefly recapitulate the basic method. Our sample is NUV 
selected in the GALEX Medium Imaging Survey (MIS). The MIS/SDSS DR4 co-sample 
occupies 524sq. deg. of the north galactic polar cap and the southern equatorial strip. Our 
sample is cut as follows: 1) NUV detection, nuv_weight > 800; 2) SDSS main galaxy sample, 
Zconf > 0.67 and specclass=2; 3) 14.5 < r < 17.6, 16 < NUV < 23.0; 4) nuv_artifact<2; 5) 
field radius less than 0.55 degrees; 6) 0.02 < z < 0.22. We generate the volume correction 
using a simple V max approach, with joint selection in N UV and r mag n itude . The NUV and 
r-band magnitudes are k-corrected to z = 0.1 using the lBlanton et al.l (120031 ) code. We have 
made no evolutionary correction. 



As shown by lWyder et al.l (120061 ) the Hess function displays a well-defined blue sequence 
and red sequence. The minimum NUV — r color between the two sequences varies slightly 
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with M r , but a good overall choice is NUV — r ~ 4. More precisely, unless otherwise stated 
we use the bin with 4.0 < NUV — r < 4.5, or NUV — r ~ 4.25. As we discuss in section 
6.1.5, the results do not depend significantly on this choice. 



3.2. Extinction Corrections 



The NUV-r color o f galaxies, and t o a le ss er extent M r , i s affe cted by dust extinction as 
well as stellar age (e.g. ] Johnson et al.l (120061 ); I Wyder et al.l (120061 ). We therefore construct 
an extinction-corrected Hess function in order to determine the number density of transition 
galaxies. We do not use the method of Johnson et al. because we use the combination 

of D„(4000) and NUV-r later as a transition rate metric. Kauffmann et al. (12003) measured 
dust extinction in the z-band (A 2 ) using the continuum fluxes. iBrinchmann et al.l (120041 ) 
measure extinction using the Balmer decrement. Using the derived A z is not ideal because 
they were derived by comparing D n (4000) , R5 \ and the continuum fluxes to a large suite of 
possible star formation histories. As we discuss in the next section, we also use these indices. 
We have used both methods to correct the Hess function and find that they yield very 
similar results (cf. $6]). We discu ss there why using the d erived A z is a reasonable approach . 
Our baseline correction uses the iKauffmann et al.l (120031 ) A z . Using a lCalzetti et al.l (119941 ) 
extinction law, we find that A r = l.QA z and Aj<ujv = 3.8A Z . 

We show evidence that the extinction correction is working reasonably in Figure HJ 
This shows the Hess function overlaid with the mean minor to major axis ratio for each 
bin. Without extinction correction (left panel), there is a clear systematic trend of lower 
inclination galaxies falling in redder NUV-r bins, since lower inclination galaxies will exhibit 
higher dust extinctions. After the correction (right panel), this trend is almost completely 
eliminated. 

As we see below, the extinction correction has two effects. It reduces the volume density 
of galaxies in the transition zone. Moreover it significantly reduces a population of transition 
color galaxies showing younger than expected spectroscopic ages. 



3.3. Star Formation History Indicators 



We use D„. (4000) and H^ as calculated and employed for the SDSS spectroscopic sam- 
ple by IKauffmann et all (120031 ) and available as the MPA/JHU DR4 Value-Added Catalog. 
H^ is corrected for nebular emission. We calculate the mean and standard deviation of these 
quantities (when measured) in each color-magnitude bin of the uncorrected Hess function. 
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In Figure [3] and [2] we show the mean values of D n (4000) and US a in each color-magnitude 
bin (shown by the bin's color), with the density contours of the Hess function superimposed. 
We show this for both the observed colors and the extinction-corrected colors. Average val- 
ues of the two indices transition smoothly from younger to older characteristic ages moving 
from the blue to the red sequence at a given M r . For a given NUV — r, ages decrease slowly 
with decrease r-band luminosity. The relative dispersion of the index within each bin is also 
shown. Typically, the dispersions are ~0.3 for D n (4000) (of the full range of about 1-2.5) 
and ~ 2 for H^ (of the full range of -2.5 to 7 in the mean). 

The dispersion in both quantities is higher in the transition zone. This is illustrated 
in Figure HI In Figure [5] we show that the variance in B.5a and D n (4000) is correlated 
as would be expected (and is not due to measurement error). We show these for both the 
uncorrected and extinction-corrected Hess function. The dispersion is somewhat lower in 
the extinction-corrected case, as some of the lower stellar age galaxies have been shifted to 
bluer NUV-r bins. The dispersion suggests that we find a range of star formation histories 
in a given NUV-r, M r color-magnitude bin, even after correction for extinction. We explore 
this issue in the next section. 



4. Quenched Star Formation Models 

We seek to explore the rate at which galaxies are transitioning from the blue to the red 
sequence. In order to complete this analysis, we need to determine the color derivative dy/dt 
in the transition region. We also hope to gather additional evidence that our star formation 
history models are a plausible description of the transition galaxies. 

Our approach is to use a very simple, parametric model. We assume that star formation 
on the blue sequence proceeds at a constant rate, and is then quenched at a time t q . We 
allow a single free parameter, the quenching rate 7. Thus the star formation history we 
assume is simply: 



M*(0)(t < t q 

M*(0) (1 



,-7(t-t<j) 



(4) 
(5) 



We use the iBruzual fc Charlotl (120031 ) stellar synthesis models to calculate the evolu- 
tion of the spectral energy distribution. We assume solar metallicity and Salpeter initial 
mass function. The metallicity assumption is reasonable because, although there is a mass- 
metallicity variation, the bulk of the mass flux comes from masses 11 < logM < 10, where 
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metallicities are near solar. The small variations around this produce even smaller variations 
in the model tracks. We also use unextinguished models. We assume that t q = 10 Gyr, which 
is long after all the indexes and colors we consider have reached stable values. We examine 
below the impact of this assumption. 

The evolution of NUV-r for the quenched models is graphed in Figure El for a selection of 
quench rates: 7 = [20,5,2,1,0.5] Gyr -1 , or r = I/7 = [0.05,0.2,1,2] Gyr. The transitional 
color NUV-r=4.25 is reached at different times after the quench time t q , increasing with I/7. 

In order to calculate the indices H5a and D n (40 0) , w e use the instantaneous stel- 



lar burst index evolution reported by iKauffmann et al.l (120031 ) and compute the continuum 



wei ghted average for each star formation history. We find that the indices calculated with 



the iPicklea ( 119981 ) stellar library are a better fit to the observed indices at later ages (partic- 
ularly B.6a ), so we use these. The evolution of 115 a and D n (4000) for the quenched models 
is shown in Figures U\ El This illustrates that for our assumed star formation history, and a 
fixed NUV-r color, R5a and D n (4000) are measures of the quench rate 7. Faster quench 
rates yield younger characterisitic stellar ages from R5a (more positive) and D n (4000) (less 
positive), at a fixed NUV-r color. 

In Figure [9] we show the evolution of vs. D n (4000) for the various quench rates, 
and the points on the locus when NUV-r=4.25. In Figure [101 we add the galaxies from the 
extinction-corrected Hess function bin 4<NUV-r<4.5 and —22 < M r < —21 (from Figure 
E]). Again, we have made no signal-to-noise cuts. It can be seen that the model indices are 
a reasonable representation of the index data and that a range of 7 can potentially account 
for much of the dispersion in a given color-magnitude bin. 

Finally, in Figure [TTl we show the color derivative dy/dt vs. color for the quenched 
models. At y — 4.25, the derivative varies from a high of ~ 3 Gyr -1 to about 0.4 Gyr -1 
for 7 = 0.5 Gyr -1 . Not suprisingly, dy/dt(4.25) monotonically tracks 7. 

We remind the reader that in the spirit of an exploratory investigation we are using very 
simple star formation history models to track the evolution in D n (4000) , H5a , and NUV-r 
space. We are also assuming solar metallicity. The star formation history assumption is 
further investigated below. We defer a more complex and realistic approach to the future. 



5. Constraints on the Mass Flux 



The machinery is now in place to calculate the transitional mass flux. For all calculations 
we use the color bin 4 < y < 4.5. We do this in three ways: 
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1. Use the uncorrected Hess function, calculate a mean value of H5 A and D n (4000) for 
each M r bin, infer the corresponding dy/dt from Figures [TU1 and [TTl and use Equation 
[3] to add up the total mass flux. 

2. Use the same methodology as 1. but use the extinction-corrected Hess function. 

3. Calculate a 7 and dy/dt for each galaxy, and an average dy/dt tot each color-magnitude 
bin. 



5.1. Method 1 and 2 



In Tables ED and we tabulate the values of 4>( M r ,y = 4.25), H5AM r ,y = 4.25) 



-D n (4000)(M r , y = 4.25), 7, dy/dt, stellar mass (from iKauffmann et al.l (120031 ) available in 
the MPIA/JHU Value-Added Catalog), and mass-flux for the uncorrected and extinction- 
corrected cases. The quench rate 7 (and the corresponding dy/dt) are found by finding the 
y=4.25 ( D n (4000) , R5 A ) model pair that comes closest to the observed ( D„(4000) , R5 A ) 
pair. Distance in the ( D n (4000) , H5 A ) plane is calculated by normalizing each index to 
its full range, viz. RS A „ = H5 A /12, and D n (4000) [norm] = D n (4000) /l.O. 

In the uncorrected case, the total mass flux is p T = 0.095 M yr" 1 Mpc" 3 (the sum 
of the right-hand column). The majority of the mass flux comes from galaxies with —22 < 
M r < —20. Lower mass galaxies have younger stellar ages, faster quench rates and larger 
dy/dt. However, their contribution to the total mass is decreased. Note that each mass flux 
entry in the table is obtained from p T = 2M(M r )<p(M r} y = A.25)dy/dt[j(D n (4:000), HS A )], 
where the factor of two accounts for the fact that the Hess function NUV-r bins are 0.5 
wide. 

In the extinction corrected case, the total mass flux is = 0.028 M Q yr _1 Mpc -3 , a 
factor of two lower. Two effects cause this: the number density is lower (since some star 
forming galaxies are moved to a bluer NUV-r bin, and the characteristic ages are somewhat 
older (for the same reason). The distribution shifts to higher M r bins, with the major 
contribution coming from —23 < M r < —20.5. This second method is likely to be applied 
in deep galaxy samples for example using stacked spectra to study the mass flux evolution. 

It is interesting to examine the distribution of galaxies vs. 7. We have binned the 
galaxies into five 7 bins, shown in Figure [121 In Figure [13] we show the color derivatives 
dy/dt, the fraction of galaxies in each bin, and the fraction weighted by the color derivative 
dy/dt which compensates for the briefer residence time of fast quench objects. While the 
unweighted distribution shows a preponderence of slow-quench galaxies, the weighted distri- 
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bution shows that the underlying distribution is fairly uniform. Galaxies with faster quench 
rates will spend less time in the transition region, and the fact that the majority of galaxies 
are found to have low quench rates in the transition zone does not mean that the majority 
of galaxies undergo slow quenches. 



5.2. Method 3 

Since there is considerable variance in the age indicators in a given color-magnitude 
bin, and the non-linear dependence of dy/dt on the various age indices, we calculate the 
7 and dy/dt for each galaxy, and then determine the average dy/dt and its error due to 
the measurement error of D„(4000) and , within a particular color-magnitude bin. 

We perform this average for each M r bin, and display the results in Table [31 The color 
derivatives and resulting mass fluxes are ~15% larger. The resulting summed mass flux is 
Pt = 0.033 ± 0.001 M yr _1 Mpc -3 . The small reported error is the formal error due to 
the Poisson error in the Hess functionand the index measurement error. The errors in the 
individual M r bins are given in Tabled and are typically ~ 10 — 25%. 



5.3. Comparison to Mass Growth of the Red Sequence 



Two recent major studies. iBell et al.l (j2004bl ) using COMBO-17, and lFaber et al.l (120051 ) 
using DEEPII and COMBO-17, have conclu ded that the red se quence has grown significantly 
over the period z=l to z=0. In particular, iFaber et al.l (120051 ) measured an increase in the 
luminosity function normalization of A<fto = 0.56 ± 0.09 dex since z = 1.0 and A0o = 
0.36 ± 0.09 dex since z = 0.8. 



We convert this to a mass flux using the following relationship: 

■ ( M \ ^10 Alog *° - 1\ 1 

Pr^Jb(RS) 1 — 1 



1O Alog0 o J t(z 



(6) 



0.1) -t(z) 

Faber et al.l (120051 ) estimate that the mass-to-light ratio = 6 for today's red-sequence 
galaxies, and the red sequence B-luminosity density at z = is js = 10 7 7 L yr" 1 Mpc~ 3 . 
This gives a mass flux of p R = 0.034(0.031) M yr" 1 Mpc~ 3 for < z < 1 (0 < z < 0.8), 
both with a minimum error of ±0.1 dex. By doing this we assume that the mass flux is 
constant over the redshift interval considered, in order to compare it to the blue sequence 
mass flux obtained "today" (where for today we use z=0.1). 



We compare these results to the derived blue sequence mass fluxes using Method 1 and 
Method 3 in Figure HH The horizontal bars indicate the redshift range over which the flux 
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has been averaged. For reasons that we discuss in the next section, we consider the blue flux 
measurement as an upper limit. It is intriguing that the blue sequence and red sequence 
mass fluxes are in rough agreement. This comparison makes sense only if we are measuring 
mass fluxes at comparable epochs (in general we are not) or i f the mass flux is constant 
over time (our principal assumption). The overall conclusions of iFaber et al.l (120051 ) are not 
inconsistent with a constant red sequence mass growth rate over < z < 1, and although 
there is some evidence for a slightly increased rate for 0.8 < z < 1, the average rate obtained 
over both redshift intervals is within 10%. We return to this issue in the section EJ 



5.4. Comparison to Mass Loss from the Blue Sequence 



To calculate the blue sequence mass flux, we can use the change in the blue lu minosity 
densit y, but we also need to account for evolution in the mass-to-light ratio. Blanton iBlanton 
(120061 ) has shown by comparing SDSS and DEEP2 luminosity and color-magnitude distribu- 
tions that the blue sequence becomes redder by {u — g)o.i[z — 0.1] — (u — p) L0 [z — 1] ~ 0.3 
and luminosity decreases by ~ 1 mag. This is consist ent with an incre ase mass-to-light ratio 
from (M/L B )[z = 1) ~ 1.3 to (M/L B )[z = 0.1] = 2.5 iBell et all (l2003h . Then, with 



Pb 



Jb(BS) (g ) [z = 1]-Jb(BS) (£) [z = 0.1] 



t(z = 0.1) -t(z) 



(7) 



and j B (BS)[z = 1] = 10 8 ' 2 and j B (BS) [z = 0.1] = 10 7 ' 8 , for b lue-sequence ga laxies gives 
—p B ~ 0.01 M Q yr _1 Mpc~ 3 . This mass flux is low, as Blanton iBlantonl (120061 ) concluded, 
at least a factor of four below our measured transition flux. This conclusion is n ot changed 
by ex tinction, since the extinction vector lies along the (M/L B ) vs. u — g locus (IBell et al. 
2003h . 



How can galaxies be populating the transition zone without depleting the mass of the 
blue sequence? As suggested by the anonymous referee, there is significant mass growth on 
the blue sequence due to ongoing star formation. If we integrate the star f ormation history 
derived for example by GALEX measurements ( Schiminovich et al.lboOol ) and then (very 
crudely) divide by the age as in 



PSF 



j z z ;M*(z)dz 

t(z ) - t(zi) 



we find that p SF ~ 0.037 ± 0.006 M yr" 1 Mpc" 3 for z 
0.031 ± 0.006 M yr" 1 Mpc~ 3 for z x = 0.8. 



0.1 and z\ = 1.0, and psf 
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Remarkably, it appears that on-going star formation balances the mass flux off the 
blue seq u ence. In fact, mass continuity, the constancy of the blue sequence mass found by 



Blantonl (120061 ). and the measured star formation history strongly argue that galaxies are 



transitioning off the blue sequence at the rates derived independently in this paper. It will 
be interesting to understand whether this balance is merely a coincidence. 



6. Discussion 
6.1. Issues and Caveats 

6.1.1. Aperture Corrections 

SDSS spectroscopy is obtained with 3 arcsecond fibers, corresponding to a physical 
diameter of 3 kpc to 10 kpc for 0.05 < z < 0.2, where the bulk of our sample lies. A possible 
concern is that this leads to aperture bias. In our case, we might worry that disk light is 
underrepresented or bulge overpresented in the H.6~a and D n (4000) indices. This would lead 
to an underestimate of the quench rate and dy/dt. If this were the case, we would expect 
that the more distant galaxies in the sample would have higher (and lower D„(4000) ) 
as the disk is brought into the fiber aperture. However, no such trend is apparent in the 
data-there is no redshift dependence at all in the index values. 



6.1.2. Time Variations and Delays 

In comparing the transitional, red, and blue-sequence mass fluxes, we have glossed over 
time delays in the transition from blue to red sequence. These could impact the present 
measurement if there are variations in the mass flux rate over the cosmic timescales spanned 
by the DEEPII and COMBO-17 observations. We have determined typical quench rates to 
be ~1-1.5 Gyr -1 . With this rate, evolution from blue to our transitional color takes ~ 3 
Gyr, with a similar additional interval to reach the red sequence. Strictly speaking, we must 
compare the mass flux rate now with the blue sequence evolu tion 3 Gyr i n the past and 



the red sequence evolution 3 Gyr in the future. The results of iFaber et al.l (120051 ) suggest 
that the mass flux rate is roughly constant, although the rate may be higher in the highest 
redshift bin. It will be of some interest to measure the transitional flux at higher redshifts, 
to determine whether it is evolving with time. 
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6.1.3. Star Formation Histories 

We have assumed that all galaxies in the transitional region follow the quenched star 
formation history of equation HI Many other star formation histories could be considered 
which would produce different evolution in R5a , D n (4000) , color space. We maintain, 
however, that the quenched models we use lead to the highest possible color derivatives, 
dy/dt, in the transitional region. In this sense, since the mass flux we derive is proportional 
to dy/dt, our result represents an upper limit. 

We assumed that the quench occurs after 10 Gyr of constant star formation. This 
produces some galaxies (with the oldest D„(4000) and H.6~a indices) that exceed the age of 
the universe by several Gyrs. If the constant SFR period is reduced to 5 Gyr, the implied 
quench rate for a given index is lower, and correspondingly the total mass flux is reduced by 
about 20% from 0.026 to 0.021 M yr" 1 Mpc~ 3 (Method 3). 

We consider three other possible histories, a smooth exponential decay in star formation 
rate from formation (smooth), recent bursts, and a composite stellar population, consisting 
of an early burst and a low level of on-going star formation (composite) . 

Smooth. In this case we assume that galaxies were formed ~10 Gyr ago with a range 
of e-folding p arameters 7 t hat accounts for the range in present day colors, the original 



assumption of iTinsleyl (119681 ) . We show the time evolution of NUV-r in Figure [151 The time 
at which the transitional color NUV-r=4.25 is reached depends somewhat upon the value of 
7. We show the color evolution of H5a and D n (4000) in Figure | 



We show a summary index comparison of the smooth models, quenched models, and 
data in Figure [T71 The smooth models follow roughly the same locus as the slowest quench 
rate curve (7 = 0.5 Gyr -1 for the quenched model. The points on the locus where NUV- 
r reaches 4.25 is shown for both model sets. Smooth models reach NUV-r=4.25 at "old" 
index values, corresponding again to the index values at NUV-r=4.25 for the 7 = 0.5 Gyr -1 
quenched model. The figure also shows the index data (for galaxies in M r bins which 
contribute the bulk of the mass flux derived above, —22 < M r < —20). The preponderance 
of galaxies falls at higher B.5a and lower D n (4000) (above and to the left) than predicted 
by the smooth models at the transitional color. Many of the galaxies fall above the smooth 
B.5a - D n (4000) locus at all colors, where the quenched models provide a better match to 
both indices at the correct color. For the small fraction of the galaxies have indices and colors 
consistent with the smooth models, the corresponding color derivative dy/dt ~ 0.3 is very 
close to that for the low quench rate model with the same 7. In other words, for long decay 
times, the prior star formation history does not affect the colors and indices: smooth and 
quenched models give similar results, and our methodology still applies. Clearly, however, 
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the majority of galaxies are not well represented in the index-color diagram with smoothly 
declining models of age ~10 Gyr. 

Smooth models with significantly higher values of 7 can produce the transitional colors 
but only if galaxies continue to form at late times. These models can reproduce the younger 
spectroscopic ages as well, and actually lead to color derivatives and mass fluxes that are 
not dissimilar to those derived for quenched models. Late formation (rather than late dry 
merging) of a significant population of red sequence galaxies does not however accord with 
most observations. 

Bursts. A massive red-sequence galaxy suffering a new starburst due to a minor merger 
with a gas-rich galaxy will evolve from the red to the blue sequence, and then may evolve 
back to the red sequence after the gas supply is exhausted by star formation or feedback. 
Some fraction of galaxies in the transitional color bins may be in this state moving in either 
direction. As long as the starburst ramp up time is short, the transition from red to blue 
will be very fast, with little impact on the transition color number densities. However the 
subsequent decay is determined again by the "quench time" of the starburst. In general 
the complexity of possible burst histories makes the present simple treatment inadequate. 
It can be shown that the quench rate, H5a and D n (4000) evolution, and color derivative 
are similar to that of the quenched case studied above, but with quantitive differences of 
order 2-3. There are other sub tleties which could influence the derivation of a mass flux from 
the red sequence evolution of iFaber et al. J2005|)(such as mass-to-light ratio). Regardless, 



bursting previously red sequence galaxies produce additional numbers in the transitional 
colors as they return with considerably lower mass flux. We are secure then in interpreting 
our mass flux measurement as an upper limit. 

Composite Histories. An early burst superimposed with low-level, on-going star for- 
mation can produce a galaxy of any intermediate color if properly tuned. We show this in 
Figure [181 For example, a 10 11 M bulge can be combined with a 10 9 75 M constant star 
formation rate disk, which corresponds to ongoing star formation that has produced 3% of 
the total mass after 10 Gyr, produces a color which hovers around NUV-r~4 after about 10 
Gyr. Objects such as this fall on the lower right of the - D n (4000) diagram: at the 

transition color, H5a = —2.1 and D n (4000) = 1.96, rather old values. In other words the 
relatively low level of ongoing star formation does significantly perturb the spectral indices 
from that of an old population. The color derivative for this model at the transition color 
is low and negative, dy/dt ~ ~0-l- Certainly, there is evidence for low level star formation 



in nearby ellipticals (lYi et al.ll2005l ). The contribution of galaxies in this bin, as we showed 
in Figure [131 is relatively small since the color derivative is low anyway. Zeroing this contri- 
bution would not signficantly affect our result. Alternatively, we can simply note that these 
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galaxies have no impact on our result interpreted as an upper limit. 



6.1.4- Extinction 



We have performed an extinction correction to remove from the transition region heavily 
reddened blue-sequence galaxies. This led to about a factor of two decrease in the d erived 
mass flux. The extinctions were obtained for the z-band by iKauffmann et al.l (120031 ) using 
continuum fitting to the star formation histories. The NUV light-weighted stellar age in 
the quenched models is ~1 Gyr. It is therefore reasonable to use the z-band continuum 
extinctions in most of these transition color galaxies because the NUV ligh t is tracing stellar 



popul ations which have diffused away from star forming regions. The IKauffmann et al. 



( 120031 ) suite of star formation histories included smooth declining exponentials plus bursts. 
Their suite did not explicitly include the kind of quenched models we use, although the 
combination of bursts and declining models could mimic them. We examined the sensitivity 
of optical colors to the quench rate and found that for example at a fixed (transitional) g-i 
color there is very little sensitivity to the quench rate in the D n (4000) vs. plane. In 

other words the optical continuum given by smooth exponential models and that produced 
by our quenched models do not yield significantly different values of US a and D n (4000) since 
all three are tracing similar timescales. This means that A z derived from the indices and the 
optical continuum should be reasonably accurate even for quenched models. A self-consistent 
approach is certainly warrented in a future study. A small fraction of the A z extinctions are 
negative-for these we assumed zero extinction. 



As a check we repeated the calculation using extinctions obtained by iBrinchmann et a 



(120041 ) from the Balmer decrement (and additional lines), again with a lCalzetti et al.l (119941 ) 



extinction law. A significant fraction of the galaxies in the transitional color bins do not have 
lines strong enough to support this technique. For those we continued to use the continuum 
based extinctions. The results of Method 3 are virtually identica l as are the values of the 
quenching rates. This is comforting, since the lCalzetti et al.l (119941 ) law assumes a particular 
relationship between Ha and UV extinction which was derived for starburst galaxies. 

Removal of more extincted, presumably star forming galaxies would likely result in lower 
values of 4>(M r , NUV — r) and dy/dt, yielding even lower values of the transitional mass flux. 
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6.1.5. Color Variation 

We chose a transitional color bin of 4.0 < NUV — r < 4.5. Is our result dependent on 
this choice? If our approach is sensible, continuity requires that the mass flux calculated 
from neighboring color bins should be similar, as long as we are probing truly transitional 
objects. Differences would arise for two reasons. First, redder transitional colors probe 
slightly older timescales. If there has been any evolution in the mass flux, then this would 
produce a trend. It is likely that the m ass flux was larger in the past (preliminary analysis 



using AEGIS data iMartin et al.l (120061 ) indicates some evolution). Secondly, obj ects that 



reach the red sequence could experience episodic residual star formation events lYi et al. 



( 120051 ) that move them back into the transition region or the bluest part of the red sequence, 



increasing the apparent mass flux. 

We have repeated the analysis, using Method 3, on two neighboring color bins. The 
results for NUV -r =3.75, 4.25 [baseline], and 4.75 are, respectively p T = 0.026, 0.033, 0.038 
M yr _1 Mpc -3 . The trend is consistent with both of the effects discussed above. However, 
the trend is modest, suggesting that to first order the basic method is sound. 



6.2. Transition Region Galaxies, AGN Feedback, and Mergers 

What kind of galaxies occupy the transitional color region of the Hess function? We 
show in Figure [T§1 a small atlas of transitional galaxy thumbnails from the SDSS. The galaxies 
are organized by mass and H5a • 

One interesting and suggestive property is the prepo nderence of AGNs. W e have used 



the value added dataset of AGN properties generated by iHeckman et al.l (120041 ) to explore 
this. 



We consider two ensemble properties. IHeckman et al.l (120041 ) calculated the luminosity 



in the OIII] line for galaxies that showed line ratios consistent with AGN or composite 
AGN/star formation ionizations. In Figure [201 we show the number fraction of galaxies with 
some AGN emission in each color- magnitude bin of the Hess function. In Figure |2"T1 we show 
the mean value of L(OIII}) for the AGNs in each bin. In each of these cases, although it 
does not make a signficant qualitative difference, we have used the extinction-corrected Hess 
function. We assume that AGN continuum in these Seyfert II objects is neglible even in the 
NUV. 



1 If this is not the case, some galaxies at transitional colors are red sequence "interlopers" . In principal 
the mass flux could be even lower when these objects are removed, although the correction could also add 
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Figure [20] is particularly striking, showing that the AGN fraction rises precisely in the 
transitional color region between the blue and re d sequence. This is of course extrem ely 
suggestive of an important role for AGN feedback (jSpringel. DiMatteo fc Hernquistll2005l ) in 
quenching star formation more rapidly than simple gas exhaustion. 

In the context of this work, we can ask whether the derived model parameters, in 
particular the quench rate, is correlated with AGN properties, such as luminosity. In Figure 
I2"2"l we show for each galaxy in the range —23 < M r < —20 L(OIII}) vs. B.5a and the derived 
quench rate 7. There is considerable scatter. Linear regression gives for log L(OI 1 1]) = 
alog7 + b, a = 0.37 ± 0.05, with a correlation coefficient of only 0.32. The trend is most 
evident in Figure [23], where the median luminosity is calculated in five 7 bins. 



There is a modest trend of LiOIII]) increasing with 7. Luminosity nominally traces 
accretion rate, which is likely to be higher in the gassier galaxies nearer the blue sequence. 
Higher H5a galaxies have more recently lost their gas fueling star formation, and there may 
be a time delay for that to starve the central engine. Thus the strongest statement that can 
be made is that galaxies with AGNs of higher luminosity tend to show higher quench rates. 



Heckman et al.l (12004 ) computed a specific growth rate for central black holes by using 



LiOIII]) to obtain a rough bolometric lu minosity and the central velocity dispersion and 
the M — a relation (ITremaine et al.ll2002l ) to obtain the black hole mass. The result is 



M 



BH 



M 



Ibh = 10 w - 8 L(Om])a 200 4 Gyr 



-1 



(9) 



BH 



which we compared with the derived quench rate in Figure [2H There is a clear trend 
of increasing 'Jbh with increasing 7. But there is considerable dispersion, and the current 
black hole growth rate is typically much less than the quench rate. Most of the AGNs i n 
the transition region are accreting far below the Eddington limit (IHeckman et al.l 12004 ). 
This suggests that if the processes of black hole growth and star formation are linked, the 
timescale for the quenching of black hole growth is considerably shorter than the quenching 
of star formation. 



Springel. DiMatteo fc Hernquistl (120051 ) modeled the effects of AGN feedback on the star 



formation history of a post-merger galaxy. Using their star formation rate histories, we have 
modelled the color and index evolution in the same fashion as our simple quenched models. 
In order to make a consistent comparison, we have assumed that their star formation history 
follows 10 Gyr of constant star formation rate at the starting rate. Figure [25] shows the 



some galaxies from bluer bins. 
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star formation history. Figure [261 shows the color derivative dy/dt for the four models (3 
with AGN feedback and various [271 shows the loci of models on the R5a - D n (4000) plane, 
compared with the quenched models and the data as in Figure [TDJ The crosses and diamond 
(for the no feedback model) indicate the points when NUV — r = 4.25. Three of the models 
have black hole feedback (and three values of dark halo mass), while the fourth shows the 
most massive object without the feedback effects. Interestingly, there is reasonable overlap 
between the numerical models and the simple quenched model. These models exhibit a 
slightly higher value of . A partial reason is that the numerical models have a large 

starburst just prior to quenching, which tends to elevate H5a (and lower D n (4000) ) at 
a given NUV-r color for the fading population. The color derivative for these models at 
the transitional color is dy/dt ~ 3 Gyr -1 with feedback and dy/dt ~ 0.6 Gyr -1 without. A 
moderate fraction of galaxies show the combination of color, H5a and D n (4000) predicted by 
the merger /feedback models, while many more show the older indices predicted for mergers 
and no feedback. But a sizeable fraction show even older indices inconsistent with the 
numerical histories. Using these tracks to interpret our data leads to lower values of dy/dt 
for a given combination of index at our transitional color, and a correspondingly lower mass 
flux than we derived above. 

Roughly speaking, about 8% of galaxies (bin one of Figure[13]) show properties consistent 
with the rapid quenching in these major merger models with AGN feedback. Some 10- 
20% show properties like the mergers without feedback, and the balance show even slower 
quenching rates. The majority of galaxies do not exhibit the rapid quench rates predicted by 
strong AGN feedback. However, it is important to remember that rapid quenches transition 
rapidly. From Figure [TBI (last column), we see that when we correct for the residence time 
in the transition zone, 50% of galaxies have quench rates 7 > 5 Gyr -1 . 

Finally, we can compare the se statistics to m e asure ments of the cosmic merger rate. For 



the redshift range 0.4 < z < 0.8, IConselice et al.1 (120031 ) obtain a stellar mass accretion rate 



due to major mergers (Mb < —19) of 1.8 x 10 3 M Q Mpc 3 Gyr 1 , which is about 7% of 
the total mass flux we measure. 



7. Summary 

We have introduced a new quantity, the mass flux density of galaxies evolving from the 
blue sequence to the red sequence. We developed a simple technique for constraining this 
mass flux that exploits the excellent separation of red and blue sequences in the NUV-r band 
Hess function. We used the volume corrected number density in the extinction-corrected 
UV-optical color-magnitude distribution, the stellar age indexes B.5a and D n (4000) , and 
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a simple prescription for spectral evolution using a quenched star formation history. Our 
estimated mass flux, pr = 0.033 M yr" 1 Mpc" 3 , although strictly an upper limit, compares 
favorably with estimates of the average mass flux that we make based on the optical and 
UV data. Galaxies in the transition zone are preferentially AGNs, although we find at best 
weak evidence that the quench rates are higher in higher luminosity AGNs. 

We note that our technique could be applied with alternative evolutionary clocks, such 
as morphology. Our approach using a single color and spectral index was simple, but a more 
refined technique would use the entire spectral energy distribution in addition to any relevant 
spectral indices. Our simple approach is fairly easy to apply to deep galaxy surveys, future 
work will study the evolution of the mass flux as well as its dependence on galaxy density 
and environment. 
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Table 1. Mass Flux Table Method 1 



M r 


logM* 
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# 
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£>n(4000) 


7 


dy/dt 
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-23.25 


0.00 


3 


,30e-08 


1 


-1.13 


1.73 
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0.93 
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0.0000 


-22.75 


11.38 


1 


.39e-06 


42 


-0.24 


1.78 


1.25 


0.96 


5.34e-06 


0.0006 


-22.25 


11.23 


8 


,80e-06 
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0.27 


1.75 
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1.29 


4.56e-05 
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-21.75 


11.07 


3 


22c-05 
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-0.02 
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1.19 
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0.0089 
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6 
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1.58 
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0.0156 
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10.65 


9 


53c-05 
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6.76c-04 


0.0150 


-20.25 


10.47 


1 


25o-04 
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3.81 
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1.14e-03 
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-19.75 
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1 


33o-04 
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6.54 
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0.0143 
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1 
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64c-04 
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2.35e-03 


0.0069 
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1 
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1 
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6 
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0.0031 



Table 2. Mass Flux Table Method 2 
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Table 3. Mass Flux Table Method 3 
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Fig. 1. — LEFT: Hess diagram with colors indicating mean b/a (minor to major axis ratio) 
for galaxies in each color-magnitude bin. In this and subsequent Hess plots the fractional 
dispersion in the quantity displayed as a color is given in each bin. The fractional dispersion 
is the actual dispersion normalized by the range given in the color bar, i.e., 0.4 is 40% 
fractional dispersion over a range in b/a of 0.28, or 0.11 This diagram shows that extinction 
in edge-on galaxies produces a significant bias in the diagram. RIGHT: Extinction corrected 
Hess diagram (using continuum derived A z ) with mean b/a in each bin. Bias observed in 
left hand plot is now virtually eliminated. 
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Fig. 2. — Contours show the volume corrected NUV-r, M r color-magnitude diagram (Hess 
function). Contours are spaced logarithmically in 10 even steps from <fi = 10~ 5 — 1CT 3 (0.2 
dex per step). Colors indicate the mean value of US a in each color- magnitude bin. The 
numeric value in each bin is the fractional standard deviation of in each bin, or the 
standard deviation divided by the typical range of . LEFT: no extinction correction; 
RIGHT: extinction correction. 
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Fig. 3. — Contours show the volume corrected NUV-r, M r color-magnitude diagram (Hess 
function). Contours are spaced logarithmically in 10 even steps from = 1CT 5 — 1CT 3 (0.2 
dex per step). Colors indicate the mean value of D n (4000) in each color- magnitude bin. The 
numeric value in each bin is the fractional standard deviation of D n (4000) in each bin, or 
the standard deviation divided by the typical range of D„(4000) . LEFT: no extinction 
correction; RIGHT: extinction correction. 
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Fig. 4. — Plots show variation of standard deviation from color-magnitude bin means of 
E5 A (LEFT) and D„(4000) (RIGHT), vs. NUV-r. Symbols give variation with M r , with 
M r = -22.25, -21.75, -21.25, -20.75, -20.25 given by symbol: plus, cross, dot, diamond, 
and triangle. 
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Fig. 5. — Covariance of R5a and D n (4000) in a single color-magnitude bin (M r =-21, 
NUV-r=4.25). LEFT: no extinction correction, RIGHT: A z extinction correction. 
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Fig. 6. — Quenched star formation history, NUV-r vs. time after quench. Linestyle indicates 
quenching rates: solid: 7 = 20 Gyr -1 , dotted: 7 = 5 Gyr -1 , dashed: 7 = 2 Gyr -1 , dot-dash: 
7 = 1 Gyr" 1 , dot-dot-dot-dash: 7 = 0.5 Gyr" 1 , 
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Fig. 7.— Quenched star formation history, R5 A [LEFT] and D n (4000) [RIGHT] vs. time 
after quench. Linestyle indicates quenching rates: solid: 7 = 20 Gyr^ 1 , dotted: 7 = 5 Gyr -1 , 
dashed: 7 = 2 Gyr -1 , dot-dash: 7 = 1 Gyr -1 , dot- dot- dot- dash: 7 = 0.5 Gyr -1 , 
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NUV-r NUV-r 

Fig. 8.— Quenched star formation history, R5 A [LEFT] and D n (4000) [RIGHT] vs. NUV-r 
after quench. Linestyle indicates quenching rates: solid: 7 = 20 Gyr -1 , dotted: 7 = 5 Gyr -1 , 
dashed: 7 = 2 Gyr" 1 , dot-dash: 7 = 1 Gyr -1 , dot- dot- dot- dash: 7 = 0.5 Gyr -1 , 
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D„(4000) 

Fig. 9. — Quenched star formation history, H5a vs. D„(4000) after quench. Linestyle 
indicates quenching rates: solid: 7 = 20 Gyr -1 , dotted: 7 = 5 Gyr -1 , dashed: 7 = 2 Gyr -1 , 
dot-dash: 7 = 1 Gyr -1 , dot- dot- dot- dash: 7 = 0.5 Gyr -1 . Large points indicate NUV-r=4.25 
for each 7. 
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Fig. 10. — Quenched star formation history, R5a vs. D„(4000) after quench. Linestyle 
indicates quenching rates: solid: 7 = 20 Gyr -1 , dotted: 7 = 5 Gyr -1 , dashed: 7 = 2 Gyr -1 , 
dot-dash: 7 = 1 Gyr -1 , dot- dot- dot- dash: 7 = 0.5 Gyr -1 . Large points indicate NUV-r=4.25 
for each 7. Dots show data for a single color-magnitude bin (M r =-21, NUV-r=4.25). LEFT: 
no extinction correction. RIGHT: extinction corrected. 
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Fig. 11. — Color derivative dy/dt for quenched star formation history vs. NUV-r after 
quench. Linestyle indicates quenching rates: solid: 7 = 20 Gyr -1 , dotted: 7 = 5 Gyr -1 , 
dashed: 7 = 2 Gyr -1 , dot-dash: 7 = 1 Gyr -1 , dot- dot- dot- dash: 7 = 0.5 Gyr -1 , 



1.0 



1.5 2.0 

D n (4000) 



2.5 



Fig. 12. — Quenched star formation history, H5a vs. D n (4000) after quench. Linestyle 
indicates quenching rates: solid: 7 = 20 Gyr -1 , dotted: 7 = 5 Gyr -1 , dashed: 7 = 2 Gyr -1 , 
dot-dash: 7 = 1 Gyr -1 , dot- dot- dot- dash: 7 = 0.5 Gyr -1 . Large points indicate NUV-r=4.25 
for each 7. Dots show data for a single color-magnitude bin (M r =-21, NUV-r=4.25) for the 
extinction corrected case. Straight lines indicate bins used to generate Figure [T3l 
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Fig. 13. — Volume corrected fraction of galaxies in 5 7 bins shown in Figure [121 Galaxies 
have —23 < M r < —20. Solid line shows unweighted fraction in each bin. Dashed line shows 
fraction weighted by dy/dt, which corrects for the shorter residence time of fast-quenching 
galaxies. Dotted line shows the color derivative dy/dt in each 7 bin. 
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Fig. 14. — Results from this work for the transitional mass flux from the blue to red sequence 
Green dashed line shows result of method 1, green dotted line method 2, gre en solid line 
metho d 3. Red points show mass flux estimated from red sequence evolution of iFaber et al. 
(120051 ). Higher point is based on evolution over < z < 1, while lower point on < z < 0.8 
Blue point shows estimate — pb based on blue sequence evolution. 
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Fig. 15. — Smooth, exponentially decaying star formation rates yield colors NUV-r=4.25 
near 10 Gyr ages. Linestyle indicates exponential constant: solid: 7 = 0.6 Gyr -1 , dotted: 
7 = 0.55 Gyr -1 , dashed: 7 = 0.5 Gyr -1 , dot-dash: 7 = 0.45 Gyr -1 , dot-dot-dot-dash: 
7 = 0.4 Gyr -1 . 
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Fig. 16. — Smooth, exponentially decaying star formation rates B.5a and D„(4000) indices 
vs. NUV-r color. Linestyle indicates exponential constant: solid: 7 = 0.6 Gyr -1 , dotted: 
7 = 0.55 Gyr -1 , dashed: 7 = 0.5 Gyr -1 , dot-dash: 7 = 0.45 Gyr -1 , dot-dot-dot-dash: 
7 = 0.4 Gyr -1 . 
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Fig. 17. — Black lines give smooth, exponentially decaying star formation rates B.5a vs. 
D n (4000) Linestyle indicates exponential constant: solid: 7 = 0.6 Gyr -1 , dotted: 7 = 0.55 
Gyr -1 , dashed: 7 = 0.5 Gyr -1 , dot-dash: 7 = 0.45 Gyr -1 , dot- dot- dot- dash: 7 = 0.4 Gyr -1 . 
Large crosses indicate epoch when NUV-r=4.25. In red, quenched models are plotted with 
the following linestyles: solid: 7 = 20 Gyr -1 , dotted: 7 = 5 Gyr -1 , dashed: 7 = 2 Gyr -1 , 
dot-dash: 7 = 1 Gyr -1 , dot-dot-dot-dash: 7 = 0.5 Gyr -1 . Large points indicate epoch when 
NUV-r=4.25 for each 7. Grey dots are data for -22 < M r < -20. 
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Fig. 18. — NUV-r color vs. time since formation for composite stellar population model with 
final bulge-to-disk ratio in mass after 20 Gyr given by the legend. The composite model 
combines a 10 11 M r = 0.01 Gyr exponential bulge with a constant star formation rate 
disk. The model with log B/D=-1.25 produces galaxies in the transition color NUV-r=4.25. 
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Fig. 19. — Sample atlas of SDSS 45 arcsecond square thumbnails of transition galaxies. Each 
row is a fixed M r bin, brightest to dimmest going top to bottom, Images are ordered left to 
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Fig. 20. — Hess function with AGN fraction plotted as greyscale in each bin. Fraction of 
galaxies in the bin that have AGN or AGN-starburst composite lines based on Heckman et 
al. (2004). LEFT: No extinction correction. RIGHT: Extinction corrected. Contours are 
spaced logarithmically in 10 even steps from = 10~ 5 — 10~ 3 (0.2 dex per step). 



-43- 




Fig. 21. — Hess function with AGN [OIIII] luminosity plotted as color map in each bin based 
on Heckman et al. (2004). LEFT: No extinction correction. RIGHT: Extinction corrected. 
Contours are spaced logarithmically in 10 even steps from <p = 10~ 5 — 10~ 3 (0.2 dex per 
step). Numbers in each bin are number of galaxies. 
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Fig. 22.— We show for each galaxy in the range -23 < M r < -20 L(OIII]) vs. R5 A [LEFT] 
and the derived quench rate 7 [RIGHT]. Color range is 3.75<NUV-r<4.75. 
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Fig. 23. — Median of log L[oni] vs. 7 in bins of A7 = 0.5. Color range is 3.75<NUV-r<4.75. 
Absolute magnitude range is —23 < M r < —20. 
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Fig. 24. — Black hole growth rate ^bh (see text) vs. ^quench- Color range is 3.75<NUV- 
r<4.75. Absolute magnitude range is —23 < M r < —20. 
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Fig. 25. — Star formation history of simulated major mergers with and without AGN feed- 
back from DiMatteo et al. (2005). SOLID: v(rotation)=320 km/s, with feedback. DOT- 
DASHED: v(rotation)=320 km/s, no feedback. DOTTED: v(rotation) = 160 km/s, with 
feedback. DASHED: v(rotation)=80 km/s, with feedback. 
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Fig. 26. — Color derivative evolution for the 4 DiMatteo et al. (2005) models. SOLID: 
v(rotation)=320 km/s, with feedback. DOT-DASHED: v(rotation)=320 km/s, no feedback. 
DOTTED: v(rotation) = 160 km/s, with feedback. DASHED: v(rotation)=80 km/s, with 
feedback. 
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Fig. 27. — Feedback models evolution in K8a and D n (4000) plane (black lines). SOLID: 
v(rotation)=320 km/s, with feedback. DOT-DASHED: v(rotation)=320 km/s, no feedback. 
DOTTED: v(rotation) = 160 km/s, with feedback. DASHED: v(rotation)=80 km/s, with 
feedback. Large symbols show values at transition color NUV-r=4.25. Plusses show models 
with feedback, diamond shows model 2 with no feedback. Red curves show simple quenched 
models for reference. 



